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We present measurements of branching fractions and charge asymmetries for the decays B — > rjK* , 
where K* indicates a spin 0, 1, or 2 Kir system. The data sample corresponds to 344xl0 6 BB pairs 
collected with the BABAR detector at the PEP-II asymmetric-energy e + e~ collider at SLAC. We 
measure the branching fractions (in units of 10~ 6 ): B(B° -> r 7 if*°(892)) = 16.5±1.1±0.8, B(B + 
r7_R-*+(892)) = 18.9 ± 1.8 ± 1.3, B(B° -> r)(Kn)* °) = 11.0 ± 1.6 ± 1.5, B(B + -> r)(Kit)l + ) = 
18.2±2.6±2.6, B(B° -» 77^° (1430)) = 9.6±1.8±1.1, and B(B+ -> r)K* + (1430)) = 9.1±2.7±1.4. 
We also determine the charge asymmetries for all decay modes. 

PACS numbers: 13.25.Hw, ll.30.Er 



Decays of B mesons to charmless hadronic final states 
are widely used to test the accuracy of theoretical pre- 
dictions. The decays involving 77 and 77' mesons have re- 
ceived considerable attention since early predictions were 
unable to explain the data. For decays of interest in this 
paper, there have been recent calculations from QCD fac- 
torization 0,0 and flavor SU(3) symmetry 

Charmless B decays to final states with strangeness 
are expected to be dominated by b — ► s loop ( "penguin" ) 
amplitudes. The branching fraction for the decay B — > 
r/K* is expected to be larger than most similar decays 
(though not as large as B — > rj'K) due to constructive 
interference between two penguin amplitudes . 

While the decay B — > rjK* (892) has been seen previ- 
ously [I, |(| , there have been no searches for states with 
an 77 meson accompanied by if* (1430) mesons, and no 
theoretical predictions exist for these decays. In this Let- 
ter we present measurements of branching fractions and 
charge asymmetries for the decays B° — > r]K*° (892) [?J, 
B+ -> f]K*+{892), B° -> r](Kw)f, B+ -> r)(Kir)l + , 
B° -> r]Kf (1430), and B+ -> r]K* + (1430), where we 
denote by (Kw)q the + component of the Kn spec- 
trum. The charge asymmetry is defined as A c h = 
(T~ — r + )/(r _ +r + ), where the superscript on the width 
r corresponds to the sign of the i? meson or the sign of 
the charged kaon for B° decays. 

The results presented here are obtained from data 
collected with the BABAR detector (described in de- 
tail elsewhere |8() at the PEP-II asymmetric e + e _ col- 
lider located at the Stanford Linear Accelerator Center. 
The analysis uses an integrated luminosity of 313 fb _1 , 
corresponding to 344 xlO 6 BB pairs, recorded at the 
T(4S) resonance (center-of-mass (CM) energy yfs = 
10.58 GeV), and follows closely the technique described 
in detail in Ref. 0. The sample is 3.9 times larger than 
that of Ref. @. 

The K* mesons are reconstructed from K + ir° 
(K*+ +v0 ), K° s n+ (K* K + w+ ), or K+n~ (K$ +r _) final 
states. All tracks from resonance candidates are required 
to have charged particle identification (PID) consistent 



with kaons or pions. We select 77, K® and 7r° candi- 
dates from the decays 77 — > 77 (?7 77 ), 77 — > 7r + 7r~7r° (17371-), 
K® — > 7r + 7T _ and ir° — > 77. We impose the following 
requirements on the invariant masses (in MeV) of the 
particle candidate final states: 490 < m 77 < 600 for 77 77 , 
520 < < 570 for 773^, 486 < ra^ < 510 for K° and 

120 < 77i 77 < 150 for 7T°. For candidates we require 
at least 3 standard-deviation (a) three-dimensional sep- 
aration between the decay vertex and the e + e~ collision 
point. Requirements are loose for the variables used in 
the maximum likelihood (ML) fit described below. For 
the Kit system, we define a low-mass region (LMR) by 
755 < n%Kix < 1035 MeV and a high-mass region (HMR) 
by 1035 < rnx-K < 1535 MeV 0. For the K* we use the 
helicity frame, defined as the K* rest frame with po- 
lar axis opposite to the direction of the B. We define 
H = cos Oh, where the decay angle 0h is the polar angle 
of the kaon momentum in the helicity frame. For the 
LMR, we require -0.95 < H < 1.0 for K*° and K*+^ + , 

and -0.7 < H < 1.0 for K*+ + ^ . For the HMR, we re- 
quire —0.5 < TL < 1.0 for all modes in order to remove 
the region in Ti. having very large backgrounds. 

A B-meson candidate is characterized kinematically 
by the energy-substituted mass 771es = (js — p B )? and 
energy difference AE = Eb — \yfs, where (Eb,Pb) is 
the B-meson 4-momentum vector, and all values are ex- 
pressed in the T(4S) frame. Sig nal events peak at zero 
for AE, and at the B mass [lfj for ttt-es, with a resolu- 
tion for AE (toes) of 30-45 MeV (3.0 MeV). We require 
|AB| < 0.2 GeV and 5.25 < 7n ES < 5.29 GeV. 

The angle 0t between the thrust axis of the B can- 
didate in the T(4S) frame and that of the rest of the 
charged tracks and neutral clusters in the event is used 
to reject the dominant continuum e + e~ — > qq (q = 
u, d, s, c) background events. The distribution of | cos 0t\ 
is sharply peaked near 1.0 for combinations drawn from 
jet-like qq pairs, and nearly uniform for the almost 
isotropic B-meson decays; we require | cos 9t\ < 0.9. Fur- 
ther discrimination from continuum in the ML fit is ob- 
tained from energy flow in the event via a Fisher discrim- 
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inant J- that is described in detail elsewhere [6j . 

For the modes with rj — ► 77, we reject B — * if *7 back- 
ground with the requirement | cos 9j cc \ < 0.86, where 8j cc 
is the 77 decay angle defined, in the 77 rest frame, as the 
angle between one of the photons and the B direction. 

When there are multiple candidates (less than 30% of 
events 0), we choose the candidate with a value of the 
reconstructed rj mass closest to the PDG mass [loj . 

We use Monte Carlo (MC) simulations for the few 
charmless BB background decays that survive the can- 
didate selection and have characteristics similar to the 
signal. We find these contributions to be negligible for 
all modes with an 77 — > ir + Tr~ir° decay except rj^K* . 
For all other modes, we include a component in the ML 
fit to account for them. 

We obtain yields and A c h for each decay chain from an 
extended unbinncd maximum likelihood fit with the fol- 
lowing input observables: AE, tties, J 7 , rn rcs (the masses 
of the rj and if* candidates), and Ti. For each event 
i and hypothesis j (signal, continuum background, BB 
background), we define the probability density function 
(PDF), with resulting likelihood C: 

V) = VMBS^PjiAE^V^VjimiJVjm (1) 



A' 



C = exp(-^)II 



(2) 



where Yj is the yield of events of hypothesis j, and N is 
the number of events in the sample. The free parameters 
of the fit are the signal and background yields, between 9 
and 11 qq background PDF parameters (see below), and 
the signal and qq background charge asymmetries. 
We determine the contributions from if* (892), (if7r)*>, 

and #1(1430) b y fits in thc LMR and HMR - Thc nt 
in the LMR includes if* (892) and (Ktt)q signal compo- 
nents (if* (1430) is negligible in this region), with the 
fixed (K 7r)g yield determined from the result of the fit to 
the HMR. For the fit in the HMR, all three components 
are included; the if* (892) yield is fixed from the result 
of the fit in the LMR, while the (Kir)* and if 2 * (1430) 
branching fractions are free in a simultaneous fit over 
the two (four) sub-decay modes for if* (if * + ). For the 
generated (if tt)*, spectrum, we use the LASS parameter- 
ization ^3 which consists of the if,* (1430) resonance to- 
gether with an effective-range non-resonant component. 
The if2(1430) is generated as a rclativistic Breit-Wigner 
shape with known mass and width |l0j . 

For the signal and BB background components we de- 
termine the PDF parameters from MC. For background 
from continuum and non-peaking combinations from B 
decays, we obtain the PDF from (ibes, AE) sideband 
data for each decay, before applying the fit to data in 
the signal region; we refine this PDF by letting all pa- 
rameters vary in the final fit. We parameterize each 
of the functions T^ig (toes ) , V s i g (AE), Vj(J-) and the 



peaking components of Vj(m TCS ) with either a Gaussian, 
the sum of two Gaussians or an asymmetric Gaussian 
function as required to describe the distribution. For 
V s ig(H.) we use a low order polynomial. Slowly varying 
distributions (all masses, AE and Tt for continuum back- 
ground) are represented by one or a combination of lin- 
ear, quadratic and phase-space motivated functions [|J. 
The fitted qq background PDF parameters are found to 
be in close agreement with the initial values. Control 
samples with topologies similar to our signal modes (e.g. 
B — > D(Ktttt)tt) are used to calibrate the simulated res- 
olutions evaluated from MC |(J. 

Before applying the fitting procedure to the data we 
subject it to several tests. In particular, we evaluate pos- 
sible biases in the yields from our neglect of small residual 
correlations among discriminating variables in the signal 
and charmless BB background PDFs. The bias is deter- 
mined by fitting ensembles of simulated qq experiments 
generated from the PDFs into which we have embed- 
ded the expected number of signal and BB background 
events, randomly extracted from the fully simulated MC 
samples. The small biases are listed in Table [I] We mea- 
sure the correlations in the data and find them to be 
negligibly small. 

We compute the branching fraction for each decay by 
subtracting the fit bias from the measured yield, and di- 
viding the result by the efficiency and the number of pro- 
duced BB pairs. We assume equal decay rates for the 
T(AS) to B+B- and B°B°. In Table [Q we show for each 
decay mode the measured branching fraction together 
with the event yield Ys, efficiency e, and A c h- The sig- 
nificance is taken as the square root of the difference be- 
tween the value of —2 In £ (with systematic uncertainties 
included) for zero signal and the value at its minimum. 

For the LMR the measurements for separate daugh- 
ter decays are combined by adding the values of — 21n£ 
as functions of the branching fractions, taking account 
of the correlated and uncorrelated systematic uncertain- 
ties described below. 

In Fig. ^ we show projections onto m^s of subsam- 
ples enriched with a threshold requirement on the signal 
likelihood (computed without the variable plotted) that 
optimizes the sensitivity. There are substantial signals in 
all four samples. For the HMR, separation of the (Ktt)q 
and if* (1430) signals is afforded mainly by the if7T mass 
and helicity shapes; projections of these distributions are 
shown in Fig. [3 The statistical correlations between the 
two signals are —0.42 in the HMR fits to both the B° 
and B + decays. 

The largest systematic uncertainties are due to the sig- 
nal and BB PDF modeling, the fit bias correction, the 
modeling of the if it mass distribution, the neutral selec- 
tion efficiency, and neglect of interference between signal 
components. The PDF modeling error is largely included 
in the statistical uncertainty since all background param- 
eters are free in the fit. The uncertainties in the signal 
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TABLE I: Fitted signal yield Y$ in events (ev.), measured bias (see text), detection efficiency e, daughter branching fraction 
product (n^i)i significance 5 (with systematic uncertainties included), measured branching fraction B, and signal charge 
asymmetry A c h for each mode. The first uncertainty is statistical and the second systematic. 



Mode 


Y<? fev 1 




e (%) Y[Bi (%) S (a) 


B <W~ 6 ) 




^7#K +7r -(892) 


407 ± 29 


±15 


24 


26 


17.6 


18.2 ± 1.4 


0.24 ± 0.07 


riairK™^- (892) 


111 ± 16 


±13 


16 


15 


6.3 


10.9 ± 2.0 


0.12 ± 0.14 


B —y r)K (892) 










18.8 


16.5 ± 1.1 ± 0.8 


0.21 ± 0.06 ± 0.02 




99 ± 16 


± 7 


11 


13 


6.9 


18.0 ± 3.2 


0.19 ± 0.16 




56 ± 11 


± 4 


8 


8 


6.1 


25.4 ± 5.5 


-0.05 ± 0.20 


m*K*J gn+ (892) 


149 ± 19 


±12 


22 


9 


8.6 


20.5 ± 2.9 


-0.03 ± 0.13 


36 ± 10 


± 5 


15 


5 


3.8 


11.9 ±3.9 


-0.23 ±0.28 


B + riK*+ (892) 










1 Q A 


18.9 ± 1.8 ± 1.3 


0.01 ± 0.08 ± 0.02 


i y y u v / 


163 ± 25 


±17 


15 


20 


r o 
5.3 


10.8 ± 1.9 


0.14 ±0.15 


^^"(Jf+TT") 


69 ± 17 


± 9 


i n 


10 


6.0 


11.4 ±3.2 


-0.18 ±0.25 


B° -> 77(^77)5° 










o. I 


11.0 ± 1.6 ± 1.5 


0.06 ± 0.13 ± 0.02 




93 ±20 


± 9 


i n 

1U 


10 


1 'i 


19.2 ±4.5 


-0.05 ±0.21 




39 ± 12 


± 6 


7 


o 
o 


1 

o.4 


18.0 ±6.3 


0.03 ±0.29 


r)- ll K* + (Kl-K + ) 


55 ±16 


± 5 


12 


9 


3.0 


13.3 ±4.2 


0.13 ±0.25 


V3,k;+(k s tv+) 


49 ± 11 


± 3 


9 


5 


4.4 


28.1 ±6.7 


0.18 ±0.22 


B+ -> 7,(^)2+ 










5.9 


18.2 ± 2.6 ± 2.6 


0.05 ± 0.13 ± 0.02 


^^(if+O 


72 ± 17 


- 1 


18 


14 


4.7 


8.4 ± 1.9 


-0.20 ±0.23 


?73^2* (^ + 7r-) 


40 ± 13 


- 1 


12 


8 


3.4 


12.5 ±4.1 


0.23 ±0.31 


B° -> r7^°(1430) 










5.3 


9.6 ± 1.8 ± 1.1 


-0.07 ± 0.19 ± 0.02 


^^^(Tf+Tr ) 


26 ± 12 


- 1 


13 


7 


2.3 


9.1 ±4.0 


-0.16 ±0.41 


^K* + {K+-K a ) 


20 ± 8 


- 1 


9 


4 


2.6 


17.8 ±7.2 


-0.82 ±0.47 


nil K;+(K^+) 


12 ± 10 


- 1 


13 


5 


1.8 


6.4 ±4.7 


0.05 ±0.58 


V3 ^K*+(K° s tv+) 


2± 5 


± 1 


10 


3 


0.2 


0.9 ±5.1 


-1.00 ± 1.56 


B+ — T7AT* + (1430) 










3.5 


9.1 ± 2.7 ± 1.4 


-0.45 ± 0.30 ± 0.02 




m ES (GeV) m ES (GeV) 

FIG. 1: B-candidate tues projections obtained with a cut 
on the signal likelihood (see text) for (a) B° -> r)K*° (892), 
(b) B° -> v(Ktv)o° (long-dashed, red) plus B° -> ?7^:| (1430) 
(short-dashed, blue), (c) B + -> -qK* + (892), and (d) B + -> 
■n(KTv)* + (long-dashed, red) plus 5° -> -qKf (1430) (short- 
dashed, blue). Points with uncertainties represent the data, 
solid curves the full fit functions, and dotted curves the full 
background functions. 




m K)l (GeV) cos(0 H ) 

FIG. 2: Projection of the signals in the HMR, obtained with 
a cut on the signal likelihood (see text): Kty mass for (a) B° , 
and (c) B + channels; 7i for (b) B° , and (d) B + channels. 
Points with uncertainties represent the data, solid curves the 
full fit functions, dotted curves the (Ktt)q portion, dot-dashed 
curves the (1430) portion, and dashed curves the full back- 
ground functions. 



PDF parameters are estimated from the consistency of 
fits to MC and data in control samples with similar final 



states. Varying the signal PDF parameters within these 
errors, we estimate the mode- dependent uncertainties to 
be 1-4 events. The uncertainty in the fit bias correction 
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is taken to be half of the correction. We estimate the un- 
certainty from modeling the BB backgrounds to be less 
than 1 event. 

Uncertainties in the reconstruction efficiency, found 
from auxiliary studies of inclusive control samples 
are 0.4% per track, 3.0% per 77/71-°, and 1.9% for a K° s . 
Our estimate of the systematic uncertainty for the num- 
ber of BB pairs is 1.1%. Published data provide 
the uncertainties for the -B-daughter product branching 
fractions (1-2%). The uncertainty due to the efficiency 
of the cos#t requirement is 0.5%. The systematic un- 
certainty for A c h is estimated to be 2%, dominated by 
tracking and PID systematic effects [l3j. 

Since our model does not account for interference 
among the components, we assign systematic uncertain- 
ties based on the m(iC7r)-dependence of the complex 
phases measured in Ref. 12], with allowance for unknown 
process-dependent overall phases. The effect is small for 
the LMR and about 10% for the HMR. For the HMR, the 
systematic uncertainties are applied after the combined 
fit, taking sub-mode errors as correlated. 

In summary, we have presented improved measure- 
ments of the branching fractions for the decays B° — > 
r]K*° (892) and B+ -> r/K*+ (892), as well as measure- 
ments of the decays B° -> r)(Kir)f, B+ -> r)(Kir)* + , 
B° -> ?7^ (1430), and B+ -> r]K* + (U30), which had 
not been seen previously. The first two supersede previ- 
ous BABAR measurements and agree with earlier re- 
sults and theoretical predictions 

0,11 0- We also cal- 
culate the branching fraction for the resonant decays to 
ijKq (1430) using the composition of (Ktt)q from 0. We 
find B(B° -> ?7A'* (1430)) = (7.8 ± 1.1±0.6±0.9) x 10~ 6 
and B(B+ -> 77^:* + (1430)) = (12.9 ± 1.8 ± 1.1 ± 1.4) x 
10~ 6 , where the third errors arise from the uncertainties 
on the branching fraction i^g(1430) — > Kir 0] and the 
resonant fraction of (Ktt)q. 

There are no theoretical predictions for the decays in- 
volving spin-0 or 2 mesons. The measured values of A c h 
are mostly consistent with zero within their uncertain- 
ties; the value for B° — ► r/K* (892) shows evidence for 
direct CP violation. 
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